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Abstract 

<N 

J> ■ Recently, a class of theories of massive gravity has been shown to be 

ghost-free. We study the spherically symmetric solutions in the bi- 
gravity formulation of such theories. In general, the solutions admit 
■^J- ■ both a Lorentz invariant and a Lorentz breaking asymptotically flat 

<^ ■ behaviour and also fall in two branches. In the first branch, all so- 

lutions can be found analitycally and are Schwarzschild-like, with no 
modification as is found for other classes of theories. In the second 
branch, exact solutions are hard to find, and relying on perturbation 
theory, Yukawa-like modifications of the static potential are found. 
^ ' The general structure of the solutions suggests that the bigravity for- 

■ mulation of massive gravity is crucial and more than a tool. 



1 Introduction 

Recently, there has been a renewed interest in the search of a modified theory 
of gravity at large distances through a massive deformation of GR (see for 
a recent review pQ), trying to extend at the nonlinear level [2] the semi- 
nal work of Fierz and Pauli (FP) [3]. FP is denned at linearized level and 
is plagued by a number of diseases. In particular, the modification of the 
Newtonian potential is not continuous when the mass m 2 vanishes, giving a 



large correction (25%) to the light deflection from the sun that is experimen- 
tally excluded [I]. A possible way to circumvent the physical consequences 
of the discontinuity was proposed in [5] ; the idea is that the linearized ap- 
proximation breaks down near a massive object like the sun and an improved 
perturbative expansion must be used that leads to a continuous zero mass 
limit. In addition, FP is problematic as a n effective theory. Regarding FP 
as a gauge theory where the gauge symmetry is broken by a explicit mass 
term m, one would expect a cutoff A 2 ~ mg~ x = (mMpi) 1 / 2 , however the 
real cutoff is A5 = [m^Mpi] 1 ^ or A3 = (m^Mpi) 1 ' 3 , much lower than A2 [6]. 
A would-be Goldstone mode is responsible for the extreme UV sensitivity 
of the FP theory, that becomes totally unreliable in the absence of proper 
UV completion. Recently it was shown that there exists a non linear com- 
pletion of the FP theory [7] that is free of ghosts, avoiding the presence of 
the Boulware-Deser instability [5]. Then the propagation of only five degrees 
of freedom and the absence of instabilities was generalized in [9]; this was 
shown also in the Stuckelberg language in [TO] . 

Quite naturally massive gravity leads to bigravity. Indeed, any mas- 
sive deformation is obtained by adding to the Einstein-Hilbert action a non- 
derivative self-coupling of the metric g, requires the introduction of an ad- 
ditional metric g. This auxiliary metric may be a fixed external field, or be 
a dynamical one. When g is non dynamical we are dealing with aether-like 
theories; on the other hand if it is dynamical we enter in the realm of bigrav- 
ity [H] that was originally introduced by Isham, Salam and Strathdee [T2] . 
One of the benefits is that such a theories are automatically diff invariant. 
The need for a second dynamical metric also follows from rather general 
grounds. Indeed, it was shown in [13] that in the case of non singular static 
spherically symmetric geometry with the additional property that the two 
metrics are diagonal in the same coordinate patch, a Killing horizon for g 
must also be a Killing horizon for g. Thus, it seems that in order that the 
Vaishtein mechanism is effective and GR is recovered in the near horizon 
region of a black hole, g has to be dynamical. The bigravity setup naturally 
leads one to explore the possible Lorentz breaking in the gravitational sector, 
due to different coexisting backgrounds. 

After a brief discussion in section [2] on how massive gravity can be casted 
in a class of bigravity theories, in section [3] we present a detailed analysis of 
both the flat Lorentz invariant and Lorentz breaking phases for the ghost- 
free potential recently found. In section H] we study the spherically symmetric 
type I solutions where the second metric is non diagonal, and type II solutions 
with codiagonal metrics. In section [5] we compare the bigravity solutions to 
those found in Stuckelberg approach. Section [6] contains our conclusions. 
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2 Massive Gravity as Bigravity 



Any modification of GR turning a massless graviton into a massive one calls 
for additional DoF. An elegant way to provide them is to work with the 
extra tensor g^ u . When coupled to the standard metric g^, it allows to build 
non-trivial diff-invariant operators that lead to mass terms, when expanded 
around a background. Consider the action 

S = J d\^~g~ k Ml TZ+^ [M 2 vl (TZ - 2m 2 V) + L matt ] , (1) 

where TZ and TZ are the corresponding Ricci scalars, and the interaction po- 
tential V is a scalar function of the tensor XjJ = g^ a g av . Matter is minimally 
coupled to g and it is described by L matt . The constant k controls the relative 
strength of gravitational interactions in the two sectors, while m sets the scale 
of the graviton mass. The action (pQ) brings us into the realm of bigravity 
theories, whose study started in the 1960s (see [UJ for early references). An 
action of the form (JTJ can be also viewed as the effective theories for the low 
lying Kaluza-Klein modes in brane world models [UJ. An additional matter 
sector can be minimally coupled to g, see for instance [TJ]. The massive 
deformation is encoded in the nonderivative coupling between g^ v and the 
extra tensor field g^ v . Clearly the action is invariant under diffeomorphisms, 
which transform the two fields in the same way (diagonal diffs)Q Taking 
the limit n — > oo the second metric decouples, and gets effectively frozen 
to a fixed background value and diffs are effectively broken. Depending on 
the background value of g^ v one can explore both the Lorentz-invariant (LI) 
and the Lorentz-breaking (LB) phases of massive gravity. The role played 
by g^ u is very similar to the Higgs field, its dynamical part restores gauge 
invariance and its background value determines the realization of the residual 
symmetries. 

The modified Einstein equations can be written ag^l 

^ + ^=2^7? ( 2 ) 
k££ + Q£ = 0; (3) 

where we have defined Q and Q as effective energy-momentum tensors in- 
duced by the interaction term. The only invariant tensor that can be written 
without derivatives out of g and g is X£ = g^c/au [LT|- The ghost free 

1 The gauge symmetry may be further enlarged to the full set of Diffi x Diff2 by intro- 
ducing a suitable set of Stuckelberg fields [6] . 

2 Whcn not specified, indices of tensors related with g(g) are raised/lowered with g(g). 
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potential [7] V is a special scalar function of Y£ = given by 



■i 



V = 22 a n V r 



n = 



...4 



(4) 



n=0 



where are the symmetric polynomials of Y 



V = l Vi=n, V 2 = t 2 -t 2 , I/ 3 = T 1 3 -3T 1 T 2 + 2r 3 , 
V4 = t\ — 6 7f r 2 + 8 Ti r 3 + 3 r| - 6 r 4 , 

with r n = tr(Y n ). As a result we have 



where (V')£ = dV/dY» and q = detX = det(g) / det(g) . 

3 Flat Solutions 

It is interesting to study the structure of both the Lorentz invariant (LI) and 
Lorentz breaking (LB) phases, starting from the (bi)flat solutions, which 
will be the benchmark for the asymptotic behaviour in the general case. 
Generically, equations © and (j3J) admit the solutions [TU IT6] 

g = rj = diag(-l, 1, 1, 1) , g = fj = w 2 diag(-c 2 , 1, 1, 1) , (8) 

where c and u are parameters to be determined by V, when one imposes 
Q — Q — 0. Let us discuss the LI (c = 1) and LB (c 7^ 1) cases, where these 
conditions give respectively two and three independent equations. 

3.1 LI Phase 

For the fiat LI background (c = 1) the conditions are as follows: 



Having two equations for a single parameter u, this means that in the LI 
phase one fine tuning is needed. This corresponds to the standard tuning of 
the cosmological constant. 



m 2 [V6S- (V Y)$\ 
m 2 q~ 1/2 (V'YX, 



(6) 
(7) 



a = 12 u 2 (6 a^u 2 + 4a 3 w + a 2 ) , 
ai = — 6 u) (4 a^co 2 + 3 a 3 co + a 2 ) • 



(9) 
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Since the chosen potential (J4J is ghost-free, the quadratic mass term is 
automatically of the FP form 

L m = -m 2 g (h^PKp ~ h^hap) tTv 1 * , 




where we introduced the 2-component column vector: h fJiU = (h fJ-u ,h t j -u ) t con- 
taining the two metric perturbations g^ u = rj^ + h^ u and g^ v = rj^ + u 2 h^ u . 
The mass parameter is 

m g = —m u (12a 4 w + a^co + 02) . (11) 

At the linearized level, a massless spin two and spin two with mass m g are 
the propagating modes. 



3.2 LB Phase 

When c ^ 1 Lorentz symmetry is broken by the VEV of g. The conditions 
are three: 

a = — 24 oj 3 (3a 4 w + a 3 ) , a x = 6 a; 2 (8 a 4 a; + 3 a 3 ) , 

a 2 = -6 w (2 a 4 w + a 3 ) . 

Among these, one determines a; in terms of the a;, the remaning two are fine 
tunings, equivalent to setting to zero the effective cosmological constants for 
each of the two metrics. Moreover, we note immediately that c is not deter- 
mined by the above equations, and is thus a free parameter. The situation is 
to be compared to bigravity with generic potentials [16], where two equations 
determine u and c, and only one fine tuning is required for biflat solutions. 
Hence, this is a peculiarity of the potentials (@| considered. The fact that 
c 7^ 1 is not determined, suggests that the potential has some flat directions 
for the metric fluctuations. This peculiar behaviour is due to the fact that 
for the ghost free potential the lapse N, related to gu in ADM formalism, is 
a Lagrange multiplier. This is confirmed already at quadratic level. 

The expansion of the potential at quadratic order gives a generic LB mass 
term of the form 

Lmass = ^ (/ioo m o^oo + 2/ioi m i^oj - h\jm 2 hij + hliTaizhu - 2/4m 4 /i o) • (13) 

While for a generic potential we have the following matrix structure [16] 
m = A C- 2 PC" 2 
mi = P 

Hl2,3 = ^2,3 P 

m 4 = A 4 C~ 2 P 



1 -1 
-1 1 



C = diag(l,c) 



(14) 
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for the potential (J4]), the masses turn out to be 



A = A 4 = 

3 (15) 
A2 = A 3 = - m 2 (1 - c) u 3 (4a 4 u + a 3 ) . 

The mass of the spatial transverse traceless propagating mode (2 DoF) is 
proportional to 



3 

m 2 gLB = m 2 -(1 - c)uj 3 (Aa^u + a 3 ) . (16) 



For c < 1 its positivity requires (4 040; + 0,3) > while for c > 1, it requires 
(Aa^u + as) < 0. Note that in the limit c — > 1, the LB phase intersects the 
LI phase but the graviton is massless, i.e. m gLB — > 0. 

Because det(P)=0, together with the massive tensor mode there is always 
a massless one in the spectrum of metric perturbations. The corresponding 
phenomenology is quite rich, and was analyzed in [16] . The linearized theory 
can be interpreted as a diff-invariant realization of massive gravity, free of 
ghosts and phenomenologically viable (no vDVZ discontinuity is present) H 
The only propagating degrees of freedom at linearized level are the spatial 
transverse traceless tensor modes (2 polarizations for each metric) physically 
representing a massless and a massive graviton (gravitational waves) oscil- 
lating one in the other and with different speeds, resulting in a nontrivial 
dispersion relation. The possibly superluminal speed c 2 in the second gravi- 
tational sector does not lead to causality violations, because the new metric 
has the character of 'aether'. The physical consequence is that gravitational 
wave experiments become frame-dependent. 

Finally, from the analysis of [16] (see table I there), we see that since 
A = A 4 = 0, at linearized level a degree of freedom is not determined, and 
an extra gauge mode corresponding to a shift in h 00 — h 00 appears. This 
is an artifact of the linearized approximation, because we know [7] that at 
nonlinear level no additional gauge invariance is preserved (beyond the four 
diffs). Therefore we expect this mode to be determined at the nonlinear 
order. 



4 Spherically Symmetric Solutions in Vacuum 

In GR, the form of the exterior solution for a spherically symmetric self- 
gravitating body is dictated by the Birkoff theorem to be Schwarzschild. 
Since in our case the vacuum Einstein equations are modified by the presence 
of the tensors Q, Q we expect that spherical solutions may deviate from 
Schwarzschild. 

3 Theories of Lorentz-breaking massive gravity were analyzed also in |171 118] . 
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In a spherically coordinate system (t, r, 8, ip), the form of g and g is 
ds 2 = -J(r) dt 2 + K(r) dr 2 + r 2 dtt 2 , 

(17) 

ds 2 = -C(r) dt 2 + A(r) dr 2 + 2D(r) dt dr + B(r) dQ 2 . 

Because of D, in general we cannot bring both metrics in a diagonal form 
with a coordinate transformation. 

Solutions fall into two classes [12], [19]: type I for D ^ and type II for 
D — 0. Since E 1 ^ is diagonal by the choice of the first metric, then also 
(V'Y)% must be diagonal because of (j5J). The only possible source of off- 
diagonal term in RHS of would be (V'Y)%, and as a result also E v must 
be diagonal. In general, the off-diagonal components of Q and Q are of the 
form 

Q r t =Ql(x D(r) (Aa 2 rB 1/2 + 6a 3 B + ai r 2 ) oc . (18) 
Thus, we have two options: either D(r) = or B(r) = uj 2 r 2 with 

6u 2 a 3 + 4(i)a 2 + a x = . (19) 

Comparing this condition with the conditions to have the flat background 
solution (Q or (1T2"]) . we realize that in the LB scenario such a condition is 
automatically realized, while in the LI one an extra constraint is required^ 

12a 4 u 2 + 6a 3 u + a 2 = 0. (20) 

Thus, for D / 0, the conditions required in the LB and LI phases are the 
same. 



4.1 Type I Solutions 

For D(r) ^ the solution can be found analytically and reads 
J=i _2mi +Ai ^ w r) * B = u 2 r 2 , D 2 + AC = cV 



(21) 

C = cW[l-^^ + A 2 r 2 \, A 



2K~ 1 m 2 2 \ A (c 2 + 1) u 2 J - C 



r J J 2 

where m 12 and c are integration constants, w satisfies eq. (fT9l) . and 

Ai = --m 2 (ao - 12a 3 o; 3 - 6a 2 u 2 ) , (22) 
3 

2 m 2 

A 2 = — — (I2a 4 u 2 + 6a 3 u + a 2 ) . (23) 



l Note that this corresponds to a massless LI graviton, m g = see eq. (TlTj) . 
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As a result, the geometry for both metrics is dS(AdS) Schwarzschild. It is 
remarkable that the deviation from Schwarzschild which is present in those 
exact solutions for different (quite similar) potentials [14], in the form of a 
nonanalytic term r 7 , is totally absent for this choice of potentials. 

Of course, if we require the solution to be asymptotically flat, then Ai >2 = 
and the three conditions ()12p must be satisfied, for both the LI and LB 
cases]! While in the LB phase this leaves space for a massive graviton (m g = 
m gLB 7^ 0), in the LI phase eq. fl20|) forces the mass of the spin 2 mode to 



A remarkable property of the type I solutions is that independently on 
the potential JK = 1; as a result no vDVZ discontinuity is present. The 
constants m\ and rri2 are related to the total gravitational mass m to t of the 
system. In fact a study of the total energy (in the asymptotically flat case) 
showed [20J that m tot = mi + m 2 , with interesting cancellation of the mass 
screening mechanism in the case of realistic star solutions. 

The key point in deriving the solution ff2T]) is that once B = uj 2 r 2 , then Q 
and Q behave exactly like the energy momentum of a cosmological constant; 
this fact is also interesting in view of cosmological applications. 

From the linearized LI phase one expects that a combination of a massless 
and massive spin 2 modes should mediate gravitational interactions, however 
surprisingly for type I solution only the massless field is important as it is 
evident from the 1/r behaviour of the solution (pZTj) . The absence of a Yukawa 
exponential suppression of the static gravitational potential can be explained 
by the fact that when (jT9l) is satisfied, the FP mass vanishes, see eq. (TTTj) . 
As far the LB phase is concerned, the values of the masses are such that all 
scalar modes mediate 1/r instantaneous interactions |16j^l though there are 
massive tensor modes that propagate. This is the main difference between 
the LI and LB phase. 

4.2 Type II Solutions 

For D = the equations of motion are in general very difficult to solve ana- 
lytically. In general even in this case we can have LI or LB flat asymptotics. 

In this section we focus on the LI case, the discussion on the LB case can 
be found in appendix [A] 

A class of exact solutions can be found if one takes as an ansatz the form 
of B that was obtained for the type I solutions, namely B = u 2 r 2 with uj 

a Eq. (TIT))) and eqs. © together, are equivalent to eqs. (|T^|) . 

6 Note: the solution (f2lj) is valid only if 0,3 + 4 tu a± 7^ 0. In case this quantity vanishes, 
the solution (|21[) is still valid except that A disappears from the equations of motion and 
remains undetermined. Also, in this case at linearized level all LB masses (|15p vanish. 

7 From eqs. (fl"5]) one can see that the paramater A M defined in [IB] which controls the 
deviation from 1/r at the linearized level, vanishes. 




7 



satisfying (fT9]) (see also |21j). Then one gets that the metrics are conformally 
related 

~9iiv = to 2 9nv , (24) 

with 

J = l_^ + Ar 2 , KJ = 1, B = u 2 r\ (25) 
r 

where m is an integration constant and the value of A coincides with Ai of 
type I solutions, eq. ( )22|) . Again, we have a dS(AdS) Schwarzschild solu- 
tion but, differently from type I solutions, both metrics are simultaneously 
diagonal with a single integration constant. 

If no ansatz on B is given, one is unable to solve the highly nonlinear set 
of equations. Nevertheless a weak field expansion is clearly viable: using the 
consistent background g^ u = r]^ u , g^ u = u 2 r]^ u , and setting 



J = l + eSJ, 6K = l + e5K, C = tu 2 (l + e 5C) 
A = oo 2 (l + e5A), B = u 2 (r 2 + e5B) , 

from the equations of motion at order e one finds 

r 

m 2 e x a [r + X g ) 2 mi 



(26) 



5K 



5J 



rX g 

— r 

2 m 2 e A s 2 mi 



r 

SB = m 2 e~ r/Xs ' 



H 1) (\ 2 9 + r 2 + X g r) (27) 



K u 2 r 



5C 



2mi 2m 2 e r l Xg 



r KUJ 2 r 



2m x - r ,x S\ + r ) [2X 2 g (nuj 2 + l) + KUJ 2 r 2 ] 
oA = m 2 e 1 9 — - — — — 5 - 

r kuj z X g r 6 



where A^ 1 = m g ^J2{u~ 2 k~ 1 + 1). 

The solution clearly shows the vDVZ discontinuity. In fact, in the limit 
m g — > (X g — > oo) one has 

5J + 5K= m2e ~ t ^- r) -> ^, (28) 
rX g r 

which does not vanish as it does in GR. Actually the weak field expansion 
is not even well defined in the m — > limit; indeed in this limit the per- 
turbations SB and 5 A diverge. Notice however that though g^ u is singular 
in the limit X g — > oo, the associated Riemann tensor is well defined for any 
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Xg. This suggests that the singular behaviour is due to the choice of coordi- 
nates rather than to a real singularity. We leave the detailed study of this 
problem, that is related to the Vainshtein mechanism, for a future work [25] . 
The Vainshtein mechanism in the Stuckelberg approach for the ghost free 
potential was studied in [22], for different potentials see [23] . 

5 Comparison with the Stuckelberg approach 

In this section we compare our previous solutions with the ones obtained 
with a frozen auxiliary metric. Our results partially coincide with [22] . 

Massive gravity can be also formulated taking the second metric as an 
absolute flat metric, and introducing a suitable set of "Stuckelberg" fields to 
recover diff invariance. Generically, a flat metric can be written as 

~ 9[1V = E;E b „ Vab , El = 8^ a . (29) 

The four Stuckelberg fields $ a are used to parametrize the "flat" vielbein 
and physically represents the global "flat" coordinates in which the metric 
g^u is flat. In this formulation of massive gravity the action is 

S = J d 4 x^g- [Ml (U - 2m 2 V) + L matt ] , (30) 

The potential is the same of @ and the equation of motion are just the ones 
of®. 

In the spherically symmetric case one can choose coordinates such that 
<3> a = Sf, and 



(31) 



ds 2 = -F(r) dt 2 + W(r) dr 2 + 2Z(r) dt dr + P(r) dQ 2 , 
d s* = - d t 2 + dr 2 + r 2 dn 2 . 

Such a choice is sometimes called unitary gauge. 

Type I solutions. For Z ^ we have the exact solutions 

P = r 2 a 2 (a positive root of a\0? + Aa 2 a + 6a 3 = 0) , 

F = 2 _ 2M +r 2 A A _ m 2 / 2 (a 3 a -6a 2 a- 12a 3 ) ^ 
r 3a 



FW + Z 2 = a 2 f 2 , W + F = f 2 + a 2 , 



where fo and M are integration constants. As usual when fo = a we are in 
a LI phase. The requirement to have solution of the form g^ = a 2 ^^ is 

a a 3 + 3ai« 2 + 6 (a a 2 + a 3 ) = . (33) 
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When a satisfies the algebraic equation given in (|32|) . the condition for a 
conformally flat solution corresponds to A = 0. 

As a general comment, in the bigravity formulation we showed that all 
solutions of type I are of the form ( 12T1) . In the Stuckelberg formulation the 
situation is less favorable and an ansatz is required. 



Type II solutions. Here Z — 0. In this case, again we have to rely 
on perturbation theory. Expanding around Minkowski space P = r 2 + SP, 
F = 1 + 5F, W = 1 + SW, we get 



5F= ^cie _ 



r 

2dV~^ (Xg + r) 



SW=- ^ / 9 ' i (34) 
cie"^ {X 2 g + rX g + r 2 ) 



5P 



r 

-2 _ ,^2 



where c\ is an integration constant and X g = m {a\ + 4 0,2 + 6 03). 

Let us compare the above type I and II solutions found in the Stuckelberg 
approach with the ones in bigravity; this is possible by taking the limit 
k — \ 00, which freezes the dynamics of the auxiliary metric. For type I 
solutions in bigravity we have 

J -> J , K ->K, C ^c 2 u 2 , 

A^A 00 = u 2 J~ 2 [J(c 2 + l)-c 2 } , £^£00= [cW^ 2 -^)] 172 

In order to compare to the Stuckelberg approach in the unitary gauge, 
we need to change coordinates to bring the second metric in a diagonal 
Minkowski form. Taking 

dt = dt' /cu + drD/C , p 2 = B(r) , (35) 

we have 

ds 2 = — J ncw dt' 2 + K new dp 2 + 2_D ncw dt'dp + to" 2 p 2 dVt 2 
ds 2 = -dt' 2 + dp 2 + p 2 dtt 2 , 

where 



(36) 



j _ J K ,T _ ° 2 + 1 n - (17) 

'J new n , Jinew ~r '-'new o o ' J^ncw n a ■ \ / 

cr w z <r era; 4 

This gives exactly the solution (1321) with the identifications a = u~ l and 
/o = c" 1 ^ 1 . 
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For LI type II solutions the story is different. In bigravity by taking 
the limit k — > oo in the weak field solution, we find that g^ u is not flat, as 
one can see from the direct computation of the associated Riemann tensor. 
In addition, for finite X g , there is no choice of mi and 777.2 for which g^ v is 
flat. However, the simultaneous limit X g , k — > 00 leads to a zero Riemann 
tensor for g^ u when m\ = 0. Introducing T^, the EMT for matter minimally 
coupled to g^, from the above considerations and from the linearized analysis 
of [16] it should be clear that mi is an integration constant related to the 
linear combinations of the sources that couples to the massless graviton, 
while 777,2 is associated with the combination of sources that couples to the 
massive graviton. For instance for k = u = 1, we have maximal mixing and 
mi oc T 00 + Too, m 2 oc T 00 — T o. Summarizing, for type II solutions, bigravity 
in the k — > 00 limit and the Stuckelberg approach give different results. This 
is rather important for studying the Vainshtein effect which is captured by 
the type II solutions. 

6 Conclusions 

We have studied the spherically symmetric solutions for massive gravity in 
the bigravity formulation of the theory. The interaction potential that we 
used was recently shown to be ghost free. The theory admits both Lorentz 
Invariant and Lorentz Breaking flat solutions, which can be used as asymp- 
totic backgrounds. Remarkably, with this interaction potential the amount 
of Lorentz breaking, i.e. the relative speed of light in the two backgrounds, 
appears as a free parameter, not determined by the interaction potential. 

The spherical solutions falls in two separate classes: type I solutions with 
D^O, where the second metric is non diagonal and type II solutions where 
D = and the two metrics are simultaneously diagonal in the same coor- 
dinate patch. For what concerns type I solutions, we found that are always 
Schwarzschild-like and the effect of the massive deformation is equivalent, 
on shell, to a cosmological constant. This might be interesting for cosmol- 
ogy. Type I solutions do not show any modification of the static part, in 
particular the "Newtonian" potential has the standard 1/r fall-off. This is 
to be contrasted with analogous solutions found with other choices of inter- 
action potentials [H], which show a non-analytic modification with respect 
to Schwarzschild. In the LI case, which requires a fine tuning, this can be 
physically understood from the fact that at the linearized level the graviton 
mass vanishes. On the other hand in the LB phase, while the static poten- 
tial is Newtonian, the spectrum contains a massive graviton tensor with 2 
DoF [IS]), which does not get excited in the spherically symmetric configura- 
tion. The other modes do not propagate at linear level, where an accidental 
gauge invariance appears, but are expected to propagate at nonlinear level. 
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Modified gravity effects appear in type II solutions where both metrics 
are diagonal. In this case, according to the interesting work of ref. [13] two 
static, spherically symmetric, nonsingular, and diagonal metrics in a common 
coordinate system must have the same Killing horizon. These results imply 
that the Vainshtein mechanism in this case cannot take place in a black hole 
when the second metric is frozen to be Minkowski. This simple fact shows 
that bigravity is more than a tool in formulating massive gravity. Type II 
exact solutions are very hard to find. Except for a special class where the 
two metrics are conformally related, one has to rely on perturbation theory. 
One can show that the standard weak field expansion (that is equivalent to 
a derivative expansion) cannot be trusted when the graviton mass is small. 
Moreover, the solution found, in the limit k — > oo differs fron the one found 
in the Stuckelberg approach. We leave the detailed study of these solutions 
for a future work [25] . 



A LB Type II Solutions and Perturbativity 

Let us discuss now the LB case and we limit ourself to asymptotically flat 
solutions. We remark first that Type I solutions in the LB case cannot be 
derived in a standard perturbative way. In fact, suppose we try to find SD 
pertubatively, expanding at the leading order Q r t or Q[ we have 

Q\ oc SD (4a 2 u + 6a 3 tu 2 + aj . (38) 

However in the LB case, 4a 2 uj + 603 u 2 + a\ = 0, and as a result SD can- 
not be determined. The problem can be overcomed by turning to the non- 
democratic perturbation theory discussed in [H] . The idea is that the metric 
perturbations are of order e, except for SD which is of order e 1 / 2 . This choice 
is enough to capture the nonperturbative features of the solutions as shown 
in [H]. 

Let us now discuss type II solutions. Defining the graviton mass as m g = 
rrigLB (see f lT6|) ). the differential equations of the perturbations do not close; 
precisely we have: 

. T/ 2Gm 1 2 SB 
SJ = - h m — , 

2Gm 1 

(39) 



SK 



r 



SA = u 2 
SC = 



(SBY 2Gmi 
\ r J r 
m 2 SB 2Gm lC 2 u 2 



k r r 



2 
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where the fluctuation SB is not determined at leading order in the weak field 
expansion. 
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